Time— metric equivalence and dimension change under time reparameterizations 
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We study the behavior of dynamical systems under time reparameterizations, which is important 
not only to characterize chaos in relativistic systems but also to probe the invariance of dynamical 
quantities. We first show that time transformations are locally equivalent to metric transformations, 
a result that leads to a transformation rule for all Lyapunov exponents on arbitrary Riemannian 
phase spaces. We then show that time transformations preserve the spectrum of generalized dimen- 
sions D q except for the information dimension D\, which, interestingly, transforms in a nontrivial 
way despite previous assertions of invariance. The discontinuous behavior at q = 1 can be used to 
constrain and extend the formulation of the Kaplan- Yorke conjecture. 
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Recent studies of chaos in general relativity and cos- 
mology highlighted the importance of the time parame- 
terization as an extra dimension in the characterization 
of chaotic dynamics pQ . Ever since Lyapunov exponents 
and other dynamical quantities were found to depend on 
the choice of the time parameter j2] , much effort has been 
directed towards an invariant characterization of chaos 
that would avoid the difficulties imposed by this depen- 
dence [21 IH [51 EJ 13 ■ However, at the most fundamental 
level, one could instead seek to explore the freedom in- 
troduced by time transformations in order to investigate 
the dependence of the dynamical quantities on the geo- 
metrical versus temporal properties of the orbits, which 
is an important and often elusive open problem. This 
problem can be traced to the question of how dynamical 
quantities change under time transformations. 

In this Rapid Communication, we consider the dynam- 
ical effect of spatially inhomogeneous time transforma- 
tions (i.e., time reparameterizations that depend on the 
phase-space coordinates). From the perspective of the 
rate of separation between nearby trajectories, we show 
that the time reparameterizations can be identified with 
local transformations of the phase-space metric. This 
implies that all Lyapunov exponents of a given orbit are 
scaled by a common factor and the resulting Lyapunov 
dimension is invariant under time transformations. We 
show, however, that the information dimension is gen- 
erally not invariant in non-ergodic systems, illustrating 
that the identity between the information dimension and 
the Lyapunov dimension of average Lyapunov exponents 
generally does not hold in such systems; noticeably, the 
other generalized dimensions remain invariant in spite of 
their dependence on the invariant measure, which does 
change. 

Numerous physical systems can be described as smooth 
dynamical systems of the form 



dx 

lib 
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metric g, which represents the phase space of the system. 
We focus on this general class of systems and consider 
time reparameterizations of the form 



dr = r(x)dt, 



(2) 



where r is a smooth and strictly positive integrable func- 
tion. We also assume that r and r -1 are bounded away 
from zero on the asymptotic sets of the system. These 
conditions assure that r(x ,t) = J Q r(x(t'))dt' , where 
x a = x(0), is a well-defined time parameter. In the case 
of Friedmann- Robertson- Walker cosmological models, for 
instance, the proper time T and the conformal time r\ 
are related through the relation dT — adrj, where the 
dynamical variable a is positive away from cosmological 
singularities [5]. But the reparameterization ^ is not 
limited to relativistic systems, in that it can represent 
any change of independent variable; parameter r could 
be, for example, a monotonically increasing angular co- 
ordinate. 

We first note [Hj that the time reparameterization 
changes an invariant probability measure from fi to fx r 
according to 



d\x r = -j. — dfi . 
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This change applies, in particular, to natural probability 
measures, despite the fact that the orbits remain invari- 
ant and ergodicity is preserved by time reparameteriza- 
tions [9]. Physically, this reflects the fact that the trans- 
formed system evolves at different speeds and hence with 
different residence times along the orbits [TP] . 
Next we note that the Lyapunov exponents, 



X(x ,v ) = Iimsup-log||w(a;o,*)||) 

t — >OC t 



(4) 



defined on a smooth Riemannian manifold M of certain 



may change as the time is reparameterized [H [2] . Here 
v(x ,i) is the solution of the variational equation of sys- 
tem (HI for an initial condition x a and an initial vector v 
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modeling the distance between nearby trajectories, and 
||-|| is the norm induced by the Riemannian metric. The 
time transformation generally changes the length and di- 
rection of the vectors v(x a ,t) for t > 0, which are then 
denoted by v r {x Q , t(x q , t)). However, we now show that 
an equivalent change can be induced by a transformation 
of the metric. 

Specifically, we construct a Riemannian metric g on M 
such that for any two vectors v and w a in the tangent 
space T Xo M of M at x a we have 

(v(x , t),w(x , t))~ = (v r (x , t(x d , i)),w r (x , t(x d , t))) g 

(5) 

for every t in some interval / around zero. Here v and w 
(v r and w r ) correspond to the solutions of the variational 
equation before (after) the time reparameterization, and 
(•, •) stands for the scalar product induced by the metric. 
We say that g and g satisfying Q are locally related by 
the time reparameterization r at x Q . 

To proceed we consider N= dim M linearly indepen- 



dent vectors va € T x M, 



1, and denote by 



Uj(t) = y(x ,Vj,t) the solution of the variational equa- 
tion of (JlJ with yj(0) — Vj. With respect to a local rep- 
resentation (and using the summation convention), the 
variational equation reads 



dyj k = df^_ 
alt dx l 



Vi\ k = l,...,N, 



(6) 



and does not depend on the metric. The same argument 
applies to the solutions Zj(t) = z(x ,Vj,T(x ,t)) of the 
variational equation after the time transformation, 



dzj k _ df k /r 



(7) 



with Zj(0) = Vj. 

On account of this, condition ^ can be restated as 



(y l (t),y J (t))g=(z l (t),z J (t)) i 



(8) 



for every i, j — 1,...,N and every t £ I. Represent- 
ing the metric tensor g locally by G(x a ,t) = (gij(x(t))), 
Eqs. ([8]) determine the choice of a family of matrices 
G(x a ,t) along the trajectory x(t). The corresponding 
system of N 2 equations can be written in matrix form as 



Y t G{x ,t)Y t — Zt G(x ,t)Z t , 



(9) 



where < denotes the matrix transpose, G = {gij) de- 
notes the local representation of the metric tensor g, 
and where Y t = (y 1 (t) 7 y 2 (t), . . . ,y N (t)) and Z t = 
(zi(t), Z2(t), . . . , z/v(i)) are N x AT-matrices having the 
vectors Vj(t) and Zj(t) as columns. The matrix Y(t) is 
invertible since v\,... , Ujy are linearly independent and 
the linear map v a — > v(x a ,t) is invertible. The latter 
follows from the fact that the flow map x Q — > x(t) is a 
diffeomorphism on M. Therefore, ^ leads to 




FIG. 1: Time-metric equivalence. Geodesic coordinates of 



metric (111 for three different initial points and contour lines 
at distances 1/4 and 1/2 from these points. In this illustra- 
tion, all geodesic lines have length 1/2 with respect to the 
new metric. 



which is the locally related metric that we sought to con- 
struct. 

A simple example is illustrated in Fig. [T] in which we 
time transform the linear flow dx/dt = 1, dy/dt = 
with r(x,y) = e~ v . Given x(0) = and arbitrary y(0), 
the metric g that is locally related to the 2D Euclidean 
metric is given by the matrix 



(<7y) 



x 1 



(11) 



G(x a ,t) - (Z t Y t -yG(x ,t)Z t Y t 



(10) 



In the figure we show local geodesic coordinates of the 
new metric, that is, given an initial point we draw the 
outgoing geodesies of fixed lengths with respect to g. 
Note that even this simple system exhibits interesting 
properties due to the shear introduced by the time repa- 
rameterization. A few observations are in order. 

First, the metric g defined by ( 10 1 does not depend 
on the initial choice of vectors Vj, j — 1,...,N. Be- 
cause the variational equations are linear, we can write 
Y t = $(x o ,t)V and Z t = ® r (x ,T(x ,t))Vo, where 
Vq = (vx, t>2, ■ ■ ■ , •Uiv) is the A^ x A^-matrix having the 
vectors Vj as columns and where $ and <I> r are the local 
representations of the evolution matrices. This leads to 
G(x ,t) = W(x ,tyG(x 07 t)W(x ,t), where W(x Q ,t) = 
& r (x , t(x , t))^~ 1 (x a , t). Second, the metric g depends 
smoothly on the initial conditions and can be extended 
in a neighborhood of x a ; a smooth extension can be gen- 
erated for every given smooth surface of initial points 
passing through x D transversely to the flow. Third, the 
dependence of the metric g on the time t indicates that, 
in general, the metric is single- valued only on some finite 
time interval along the trajectories. For example, the in- 
terval / is generally limited by to if 

Xq IS cl periodic point 
of period to (with respect to the time t) and the function 
r is such that &(x ,to) ^ & r (x , t(x d , to)). The interval 
/ is similarly constrained by the recurrence of orbits to 
the neighborhood in which the metric is extended |llj . 
Furthermore, the time dependence of the metric can be 
eliminated in favor of a dependence on x Q and x only. 
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Therefore, our result establishes a local equivalence be- 
tween time and metric transformations on M. 

A neat implication of this time-metric equivalence is 
that, under the time reparameterization ([2]), the Lya- 
punov exponent in Eq. Q changes exclusively due to 
the transformation of the factor 1/t. The contribution 
due to the logarithmic factor remains unchanged be- 
cause the norms induced by the metrics g and g are 
logarithmically equivalent along each orbit [IT]; that is, 
for each x Q and v Q there is a sub-exponential function 
C = C{t) such that \\v(x ,t)\\g = \\W(x ,t)v(x ,t)\\ g = 
C(t)\\v(x ,t)\\ g [T5]. Therefore, the Lyapunov expo- 
nents transform as X r (x ,v ) — X(x 07 v )/A(x ), where 
A(x ) = Hindoo t(x , t)/t. This extends the result pre- 
viously derived in Ref. [1] for Euclidean phase spaces to 
the more general case of Riemannian manifolds. 

We now turn to the transformations of fractal dimen- 
sions, which cannot be accounted for by metric changes. 
The box-counting dimension Dq is purely geometrical 
and hence does not change under time reparameteriza- 
tion. The generalized dimensions D q , however, can in 
principle change for q > given that they depend on the 
measure and the measure is transformed according to 
([3]) . To analyze this dependence, we consider a positive- 
measure set of interest S (typically an attractor) and 
define the spectrum of dimensions on 5 as 

1 1 N(£) 

D M = rlimsup- log yVSk)', q>0,q£lV) 

q-l e^O log e J 



JV(e) 



Di(/x) = limsup — *— V" /i(-Bfe) log fi{B k ), (13) 
loge ^ 

where the sum is taken over the N(e) nonzero measure 
boxes Bk of edge length e necessary to cover the set 
[13l IT4l [15] . This spectrum includes as special cases the 
information dimension (q = 1) and the correlation di- 
mension (q — 2). For consistency, the measure is always 
normalized to 1 on S (with M replaced by S in Eq. Q). 
The dimensions D q are known not to depend on smooth 
transformations of the phase space. We now consider 
their behavior under time transformations. 

The first surprise is that, contrary to what intuition 
may suggest, the dimensions defined by Eq. (12 1 are in- 



variant with respect to time transformations not only for 
<7 = but also for all q ^ 1. This follows from the fact 
that ji and fi r in Eq. (|3| are absolutely continuous with 
respect to each other, i.e., both measures define the same 
sets of nonzero measure, and that r and r -1 are bounded 
away from zero on these sets. Then there exist positive 
constants ci, ci such that c\\i(Bk) < n r (Bk) < C2^(5fc) 
for every fc and e > 0. Using this in the definition (12 1, 
we obtain 



D q (ji r ) = D q (fx) for q > 0, q ^ 1, 



(14) 



i.e., the dimensions remain unchanged despite their de- 
pendence on the measure, which generally changes. 



The second surprise is that the information di- 
exhibits a distinctive behavior and 



(13) 



mension 

may change under the same time reparameterization 
To appreciate this, we first notice that Eq. (13 1 
can be written in the continuous form -Di(^t) = 
limsup e ^Qj^-j J s \ogn(B(x,e))dfi(x), where B(x,e) is 
an open ball of radius e centered at x. We then use the 
Fatou lemma to obtain 



Di(fi) < / limsup log fj,(B(x,e)) dpk{x), 

Js e-^o loge 



(15) 



where the integrand is the pointwise dimension, which we 
denote by T>^(x). A similar inequality holds for the infi- 
mum, in which case the pointwise dimension is denoted 
by V^x). This leads to 



/ V{x)d i i{x)<D l {^)< [ V^(x)d^(x). 
Js Js 



(16) 



In the remaining part of the paper we limit the discussion 
to the case P„(x) = T>^(x) = T>^(x) almost everywhere, 
a property found in many physical systems and demon- 
strated for flows with strong hyperbolic behavior |16j . 
This assures that the equalities hold in (16 1. 



The transformed information dimension is then written 



as 



Di (lJL r ) 



r{x) 
S Is rd V 



T> li {x)dn(x), 



(17) 



where we have used the measure ^ normalized on 5* and 
the invariance of the pointwise dimension. The latter fol- 
lows from ([3jl and is stated as T>^ r {x) = T>^{x) for almost 
every x. The result in (17) indicates that D\ is in gen- 



eral noninvariant when T>^(x) is not almost everywhere 
constant. 

For example, consider a system with two ergodic 
components, Sa & n d Sb, of information dimension 
Di(/j,\Sa) > Di(ijl\Sb)- For simplicity, assume that the 
original measure is evenly split between the two sets, i.e. 
h(Sa) = /J-(Sb), and that fJ,(Bk) is the same for all the 
nonzero measure boxes B^ of each set. The informa- 
tion dimension of S — Sa U Sb is -Di(a«) = ^[Dq(Sa) + 
Dq(Sb)]- Now, imagine a time reparameterization that 
changes fi\S\ uniformly by a factor < a < 2 and fi\S2 
uniformly by a factor /3 = 2 — a. The transformed infor- 
mation dimension is Di(/j, r ) = ^Dq(Sa) + "^^-D^Sb), 
which differs from -Di(^) for any a ^ 1 [2"U] , 

In the case of q ^ 1, this change in the measure con- 



tributes an additive term to log^fcM(^fc) 9 in Eq- (12 1 



that vanishes when divided by log e in the limit of small 
e, in agreement with our prediction that the other dimen- 
sions D q are all invariant. At first sight the noninvariance 
for q — 1 may seem to violate the monotonicity of D q , 
which was previously proved to hold for < q < 1 and 
for q > 1 [13], but this intuition is misleading because 
lim 9 ^i_ D q > limq_ f i+ D q whenever T>^{x) is not con- 
stant and D\ is not invariant. In our example, when 
q < 1, the contribution ^ fe ji(Bk) q from the set with the 
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D, 



1 



FIG. 2: Dimension change. The gap between the generalized 
dimensions on the left and right sides of q = 1 is invariant 
and determines the interval of possible changes for D\ under 
time transformations. The Lyapunov dimension Dl remains 
invariant, whereas the Lyapunov dimension Dl transforms as 



largest box-counting dimension dominates and leads to 
Dq<i — Dq(Sa), just as in the case q = 0; when q > 1, 
the box-counting dimension of the other set dominates, 
and D q> i = D (Sb)- The information dimension Di([i r ) 
is thus a weighted average of the dimensions on both sides 
of the discontinuity and is in general free to vary between 
lim 9 ^i_ D q and lim g ^i + D q under time reparameteriza- 
tions, as shown in Fig. [2j 

The information dimension is guaranteed to be invari- 
ant only in special cases. The most important such case 
is when [i (and hence /j, r ) is ergodic in S. Since the 
flow map x — > x(t) is a diffeomorphism and the mea- 
sure is invariant under this map, one can verify that 
V^Xo) — V^(x(t)). Then, if fx is ergodic, T>^(x) is 
constant for almost every x, and hence Z?i(/z r ) = Di(fi). 
The general condition for Di to be invariant with respect 
to any time transformation is that V^i^x) is constant al- 
most everywhere. 

It is of interest to analyze the meaning of the nonin- 
variance of D x for the Kaplan- Yorke conjecture [2TJ [55] 
and its generalizations, which state that the information 
dimension typically equals the Lyapunov dimension. Let 
the average Lyapunov exponents be \ = J Xi(x)d,Li(x), 
where Ai(cc) > \2{x) > ••• > Xn(x) corresponds to 
the ordered set of Lyapunov exponents Q at x, and as- 



sume that this definition is applied to measures that are 
not necessarily proved to be ergodic. Based on the def- 
inition operationally used in numerical experiments, the 
Lyapunov dimension can be denned as 



K 



D L (n) = K 



|Ajc + i| 



i=i 



(18) 



where K is the largest integer such that Y2j=i — 0; 
under the condition that the r.h.s. terms in (18) are 



well defined. It follows that Dl(li) remains invariant un- 
der time reparameterizations, thus violating the equality 
D L = Di when D\ changes. In the example considered 
above, the Lyapunov dimension of S = Sa U Sb is in- 
termediate between the Lyapunov dimensions of Sa and 
Sb, indicating that Dl equals D\ for at most one value 
of a. The conjecture can be re-established, however, for 
the Lyapunov dimension defined as 



K( X ) 

Xi(x)}dfi(x), 



|Ax(<r)+i(a;)| f-[ 

where K is defined as above but now at each point 
x (with the convention that the integrand is zero for 
\i(x) < and TV for \ N (x) > 0). It follows from (fl7l 



and the Kaplan- Yorke conjecture for typical ergodic sets 
that the identity D\ = Dl is expected to hold true for 
generic systems (see Fig. [2]). 

The noninvariance of the information dimension, which 
was previously surmised to be invariant [5] , is important 
as it limits the applicability of the identity Dl = D\ 
in nonergodic systems and ergodicity is a property of- 
ten difficult to verify. The invariance of other indica- 
tors of chaos established in this paper is relevant to the 
study of a range of dynamical phenomena, including spa- 
tiotcmporal chaos, and clarifies longstanding problems in 
relativistic chaos. It shows, in particular, the observer 
invariance of the often questioned chaoticity of the mix- 
master model for the early universe [3], which was first 
recognized as a chaotic geodesic flow on a Riemannian 
manifold by Chitre in 1972 in a work that made one of 
the very first uses of the term "chaos" in dynamics [TH] . 
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